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ABSTRACT 

A model rotor was mounted horizontally in the 
settling chamber of a wind tunnel to obtain perfor- 
mance and wake structure data under low climb con- 
ditions. The immediate wake of the rotor (tp < 31°) 
was carefully surveyed using 3-component particle im- 
age veiocimetry to define the velocity and vortical 
content of the flow, and the data used in a subsequent 
study to validate a theory for the separate determi- 
nation of induced and profile drag. Measurements 
were obtained for two collective pitch angles intended 
to render a predominately induced drag state and 
another with a marked increase in profile drag. A 
majority of the azimuthally directed vorticity in the 
wake was found to be concentrated in the tip vortices. 
However, adjacent layers of inboard vorticity with op- 
posite sense were dearly present. At low collective, 
the close proximity of the tip vortex from the previ- 
ous blade caused the wake from the most recent blade 
passage to be distorted. The deficit velocity compo- 
nent that was directed along the azimuth of the rotor 
blade was never more that 15% of the rotor tip speed, 
and except for the region of the tip vortex, appeared 
to have totally disappeared from the wake left by the 
previous blade. 


NOMENCLATURE 

a speed of sound, m/s 

A area of rotor disk, m 2 

c chord, m 

Cq torque coefficient, torque/ pirti 2 R 5 

Ct thrust coefficient, thrust / pnQ 2 R 4 
FM hover induced power /actual power 

M hover tip Mach number, RQ/a 

Q rotor torque, N-m 

r radial distance from vortex center, mm 

r c vortex core radius, mm 

R radius of rotor, m 

Re hover tip Reynolds number, cRQ/v 

T rotor thrust, N 

it radial velocity component, m/s 

v azimuthal velocity component, m/s 


w axial velocity component, m/s 

Voo velocity upstream of rotor, m/s 

Vup velocity of rotor tip, SIR, m/s 

6 rotor collective pitch angle, deg 

v kinematic viscosity, m 2 /s 

p air density, N/m 3 

a rotor solidity, blade area/disk area 

xp azimuthal distance from blade, deg 

vorticity normal to r-z plane, 1/s 
angular velocity of rotor blade, rad/s 


INTRODUCTION 

The performance of a helicopter rotor blade is 
influenced by two types of drag: induced and pro- 
file. Induced drag is primarily related to the con- 
centrated vortex that trails from the tip of the ro- 
tor blade. Profile drag is related to the chordwise- 
directed viscous forces that develop on the surface of 
the blade which give rise to both positive and negative 
layers of vorticity oriented along the span. Neither 
of these drag components can exist without viscos- 
ity. Induced and profile drag forces, which combine 
to produce the torque that the helicopter must over- 
come, can be reduced by cleverly tailoring the shape 
(profile, twist, planform, etc) of the blade. To fa- 
cilitate this endeavor, it is important to distinguish 
between these two sources of drag. Since neither drag 
component can be measured directly, a theory was 
developed for separately determining the value of the 
induced and profile drag components by examining 
the velocity and vortical content of the flow in the 
wake immediately behind the rotor blade (ref. 1). 
The evidence needed is believed to be present within 
the first few chord lengths behind the blade, after 
which viscous diffusion causes the layers of spanwise 
vorticity to undergo self annihilation. The purpose 
of this experiment was to gather the data needed to 
validate this theory. This paper describes the test fa- 
cility, model, flow measurement technique, and data 
acquisition and processing. A discussion of the results 
is provided along with specific findings. 
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TEST DESCRIPTION 

Test Facility- This experiment was performed 
in the settling chamber of the 7- by 10-Foot Wind 
Tunnel at Ames Research Center (fig. 1). The max- 
imum speed that could be obtained in this section of 
the tunnel circuit was about 18 m/s. Screens that are 
located upstream of the settling chamber rendered the 
flow approaching the rotor with estimated turbulence 
levels below 1.2%. 

A model rotor was placed in the center of the set- 
tling chamber and mounted horizontally with its axis 
of rotation coincident with the centerline of the test 
section. Aside from the influence of the tunnel walls 
(approximately 15 feet away) and the contraction into 
the test section (with an area ratio of 13.3:1), the ro- 
tor wake was allowed to pass through the test section 
without interference. To simulate a hover state, a cur- 
tain was drawn across the heat-exchanger section to 
prevent the wake from recirculating around the tun- 
nel and creating an unwanted climb condition. In this 
configuration, any portion of the wake reaching this 
section would be expelled through the heat-exchanger 
doors and fresh stagnant stir would be drawn by the 
rotor as required into the tunnel circuit. To simulate 
a climb condition, the curtain was removed and the 
tunnel fan was operated in a manner similar to an 
earlier experiment (ref. 2). That is, the upstream ve- 
locity was a combination of flows generated by both 
the rotor and the tunnel fan. 

Video cameras were mounted upstream and to 
the side of the rotor to monitor critical elements dur- 
ing the test (fig. 2). It is customary to infer the ve- 
locity in the test section from the differential pressure 
between the settling chamber and the test section. 
However, this procedure does not reveal the nonuni- 
form character of the flow across the test section and, 
at low speeds, is not sufficiently accurate. Therefore, 
the velocity in the test section was measured with a 
vane anemometer (1 to 130 ft/s range) attached to 
an existing tunnel traverse. In the settling chamber, 
where the velocity was considerably lower, a hot-film 
probe was used to measure velocity (0 to 8 ft/s range) 
directly upstream of the rotor at a distance of about 
28 feet from the tip-path plane. A flow seeder was 
placed sufficiently far upstream to minimize any dis- 
turbance to the flow approaching the rotor and to al- 
low enough time for the particles to spread uniformly 
over the near- wake region of the rotor to be studied. 

Model Rotor- The model used in this experi- 
ment was a 1 /7-scale of the AH-I/OLS rotor. The 
profile of this blade was based on a simple trailing- 
edge extension to the shape of a NACA 0012 airfoil 
(fig. 3). The rotor consisted of two blades, each with 
a chord of 4.09 inches and a span of 37.7 inches mea- 
sured from the axis of rotation to the blade tip (fig. 


4). Beginning at 25% radius, the twist of each blade 
was -0.27°/inch, so that the total variation in geomet- 
ric twist along the span was 7.55°. Collective pitch 
settings were referenced to the geometric pitch angle 
at 75% radius. Therefore, a collective pitch setting of 
zero meant that the blade pitch was about 5,00° at 
the root (held fixed for locations less than R/4) and 
-2.50° at the tip. The rotor hub had a built-in blade 
precone of 1.30°. 

Information related to safety-of-flight was ob- 
tained by flap- and chord-bending gages mounted on 
the blade. The drag link used to adjust the blade 
sweep was instrumented with a strain gage that re- 
sponded to the aerodynamic drag on the rotor, and 
as such, offered a redundant measurement of chord- 
wise bending. Depending on the sweep of the blade 
and the location of the center of mass of the blade, 
this gage could also be loaded by a centrifugal force 
induced moment. Strain gages were placed on the 
pitch horns that transmitted commands to the ro- 
tor blades. Reactions to these commands resulted in 
bending loads on these elements that were propor- 
tional to the torsional moments on the blades. A 
transducer for measuring the angular tilt of the rotor 
(called the teeter pot) was mounted on top of the hub. 
The rotor operator monitored the signals from each 
of these gages. Since the center of mass of each rotor 
blade was located aft of the feathering axis (aligned 
with the quarter chord), it was necessary to sweep 
the blades forward 1.20° to minimize steady chord- 
wise bending due to centripetal force. The blades 
were strobed at 2/rev to verify that the blades were 
tracking. A blade sweep of 1.20° corresponded to a 
0.19c movement of the tip. 

The model rotor was attached to a teetering hub, 
which in turn ivas mounted on the VSB-54 force and 
moment balance. The most significant element of the 
balance in the present test was the axial gage, which 
had a capacity of 350 lbs and provided a measure 
of the rotor thrust. Rotor torque was measured us- 
ing torsion gages mounted on a flex coupling located 
between the balance and the transmission, and had a 
capacity of 54 ft-lbs. The rotor was driven by a 90 HP 
electric motor through a 4:1 transmission. The com- 
plete rotor and drive system (fig. 5) was supported 
on a tripod-like structure that was rigidly attached 
to the frame of the wind tunnel. This setup placed 
the tip-path plane of the rotor about 6 feet upstream 
of the beginning of the contraction section of the tun- 
nel, well within the protective steel- wall section of the 
tunnel. 

Performance Data- The primary variables 
needed to establish the performance of the rotor are 
given in table 3. This data was sampled at 10 Hz and 
then averaged over a minimum period of 10 seconds. 
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Some cases were averaged over an extended period 
of time, sometimes as long as 30 minutes, to resolve 
equilibrium issues. All safety-of-flight data were sam- 
pled and digitally recorded at a rate of over 100,000 
samples/s and displayed at 1-second intervals on a 
monitor. 

The central region of the flow sampled with the 
vane anemometer at intervals of about 8 inches in the 
horizontal direction and 9 inches in the vertical di- 
rection. This matrix of measurements covered about 
36% of the cross section. The matrix was extended in 
both directions in a few cases to ascertain the value of 
surveying an additional 37% of the test section area. 

Table 1: Rotor Performance Variables 

Collective pitch angle 
Rotor thrust 
Rotor torque 
Air temperature 
Settling chamber velocity 
Test section velocity 


Stereoscopic PIV Method- The velocity of 
a fluid can be inferred from the motion of discrete 
particles that are suspended in the flow. This method 
is known as particle image velocimetry (PIV). A single 
thin sheet of light is normally used to define the plane 
of interest. The particle images that are recorded at 
two different times are then cross- correlated to yield 
displacement vectors (magnitude and direction). The 
time interval between the two images must be short 
enough so that the particles remain in the light sheet, 
yet long enough so that the particle displacements are 
perceptible. It is assumed that the particles are small 
enough to accurately track the flow. 

A single camera placed normal to the sheet of 
light provides the images taht are needed to produce 
a two-dimensional array of velocity vectors. If the flow 
has a significant third component (directed normal to 
the plane of the light sheet), an error due to perspec- 
tive develops in the two-dimensional array that is zero 
in the center and increases toward the image bound- 
aries. This seemingly undesirable sensitivity to out- 
of-plane motion was later exploited to derive the third 
component of velocity (ref. 3). In fact, sensitivity to 
the out-of-plane component of velocity over the entire 
image plane can be increased by intentionally placing 
the camera at an oblique angle to the light sheet. To 
uniquely determine a three-dimensional particle dis- 
placement requires two cameras that are oblique to 
the light sheet, each offering a different perspective of 
the particle motion between the two exposures (fig. 
6). This procedure is known as Stereoscopic PIV or 
3D-PIV. 


The technique demands both critical focus and 
maximum brightness of the particle image. Conven- 
tional lens mounts create a plane of focus in object 
space that is parallel with both the image plane and 
the lens plane. When the lens plane is oblique to 
the object plane (as is the case for Stereoscopic PIV), 
particles will be increasingly defocused as their dis- 
tance from the centerline of the object plane increases. 
The conventional corrective action would be to in- 
crease the depth of field by reducing the aperture of 
the lens. Unfortunately, reducing the lens aperture 
drastically reduces both the particle image intensity 
and the image resolution. The problem is solved with 
Scheimpflug focusing (refs. 4 and 5). To keep the 
object plane in focus when the lens plane is rotated 
to an oblique orientation, the image plane (where the 
camera sensor is placed) must also be rotated such 
that all three planes (object, lens, and image) inter- 
sect along a common line. 

Arrangement of PIV Equipment- Elements 
of a basic Stereo PIV system normally consist of a 
pulsed laser, optics to spread the beam into a thin 
sheet, CCD cameras to record the images, and a cali- 
bration target. Due to numerous physical constraints 
in this test, the blade was considered to be at zero 
azimuth when the trailing edge was horizontal. In- 
terest in the near wake behind the rotor therefore re- 
quired that the plane of the light sheet be horizontal 
as well, and that it should skim the trailing edge of 
the rotor (fig. 7). In order for both cameras to be 
in a forward scatter position, two opposing, copla- 
nar light sheets were needed. To avoid interference 
with the model supports, the cameras were at differ- 
ent distances away from the light sheet (nominally 
5 feet) and at different angles. The relatively low 
power of the Nd:YAG lasers that were available (120 
i nJ) meant that two such lasers would have to be 
synchronized. A 1/rev reference signal was obtained 
by directing a small continuous-wave laser beam onto 
a small mirror attached to the hub of the rotor and 
intercepting the reflected beam with a fast-response 
photo detector. 

Camera Specifications- Images were acquired 
with 8-bit charged-coupled device (CCD), cross- 
correlation cameras having a sensor array of 2K x 2K 
pixels (nominal). Each pixel measureed 7.4 /im on a 
side. The cameras were operated in a double-exposure 
mode to acquire two non-interlaced, full-frame images 
in a single frame interval. The time interval between 
images was variable between 2 /is and 30 ms. These 
are non-standard video cameras that could be exter- 
nally triggered and driven at any frequency up to 15 
Hz in double-exposure mode. A computer interface' 
provideed control over gain, contrast, black level, and 
trigger mode. Because the cameras were located at 
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different distances from the object plane, one used a 
55-mm lens while the other used an 85-mm lens. The 
lenses were remotely translated to focus on the center- 
line of the image area and the sensor (located inside 
the camera body) was rotated about its centerline to 
satisfy the Scheimpflug condition. 

Laser and Sheet Optics- A Nd:YAG laser was 
frequency doubled to provide a beam having a wave- 
length of 532 nm, a pulse width of about 4 ns, and 
power of 120 mj. The laser beam had a diameter of 
about 5 mm and a divergence of 0.60 mrad. In this 
test the laser-light sheet was formed by first passing 
the beam through a set of spherical lenses (consisting 
of -500 mm and +500 mm). By adjusting the distance 
between these two lenses, the beam waist (approxi- 
mately 1-mm thick) could be positioned in the region 
imaged by the cameras. Two cylindrical lenses (-12.7 
mm and +25.4 mm) were used to expand the beam 
in the horizontal direction. The distance between the 
sheet forming optics and the imaged area was about 
5 feet. 

System Alignment- The lasers, sheet-forming 
optics, and cameras were mounted on a common 
structure, which could be traversed in 3 orthogonal 
directions. A 6- by 9-inch rectangle was milled into 
the surface of a flat aluminum plate and used as the 
calibration target. The calibration plate was leveled 
and the light sheets adjusted so as to evenly graze the 
surface. The lens of each camera was traversed along 
its optical axis and the CCD sensors rotated until a 
uniform focus on the calibration target was achieved. 
After recording the calibration images, the structure 
supporting the entire optical system was elevated and 
translated to match the wake region of interest. 

PIV Software- This experiment was performed 
using Provision software developed by Integrated De- 
sign Tools (IDT). This software controled both data 
acquisition and the image-processing interface. This 
dual role had the advantage of providing a fully inte- 
grated calibration procedure that allowed for acquisi- 
tion, quality assurance and data reduction. The cal- 
ibration yielded all the optical parameters required 
for the accurate reconstruction of three-component 
velocity- vector fields. 

Images were processed by first covering the re- 
gion of interest with an interrogation grid. The in- 
tersections of the vertical and horizontal grid lines 
defined the centers of each interrogation window. 
The interrogation window was the smaller region in 
the reference and delayed images that were cross- 
correlated using an optimized Fourier transform. 
When the camera calibrations were applied, the grid 
in each camera view covered the exact same area in 
the flow field. 

The size of the interrogation window was based 


on the particle density and the maximum instanta- 
neous displacement of the particles. Typically, win- 
dow sizes range from 8 to 64 pixels on a side and each 
interrogation window yields one vector. The cross- 
correlation of the interrogation window in the refer- 
ence image with that of the delayed image yields a 
correlation map. The location of the correlation peak 
in space determines the local displacement in both 
magnitude and direction. 

Provision incorporates quality-assurance tests in 
the vector calculations. The processing begins with 
a first-pass correlation between the reference image 
and delayed image to determine the maximum dis- 
placement range. For the second pass the software 
enlarges the interrogation area in the delayed image. 
The amount of the enlargement is based on the maxi- 
mum displacement range determined in the first pass. 
Since the Fourier transform requires both interroga- 
tion areas to be identical, the software “adds zeros” to 
the interrogation area of the reference image to match 
the size of the delayed image area. This technique 
maximizes the probability of correlating all the parti- 
cles found in the reference image interrogation area to 
those found in the delayed image, thereby maximizing 
the statistical accuracy. 

The software then counts the number of parti- 
cle images (each composed of several contiguous pix- 
els) in the interrogation window. Ten particle pairs 
are required to contribute to a correlation map. The 
centroid of each particle is calculated, thus yielding 
sub-pixel accuracy of their position. A second-order 
curve fit is determined from the ten displacements. 
This curve fit yields a single vector, at the precise grid 
point location, whose error is reduced by a factor of 
0.3 (1A/AT, where N = 10 particles) over the straight 
fast Fourier transform (FFT) of the same area (ref. 
6 ). 

If there are not enough particles in a given inter- 
rogation region, the program will automatically en- 
large the interrogation area of the reference image. If 
there are ten or more particles in this enlarged area, 
then the correlation proceeds in the same manner de- 
scribed above. The data for this grid point is con- 
sidered to be a recalculation. If ten particles are not 
counted, the software will expand the interrogation 
area incrementally. It will repeat this process until 
the interrogation area expands to a limit of 64x64 
pixels. If the ten-particle threshold is not met for the 
largest area, then velocity values are interpolated for 
that grid point using a second-order curve fit based on 
nearest-neighbor values. This vector is also counted 
as a recalculation. The software tracks the number of 
these recalculations and displays that number after 
each camera view is processed. If the number of re- 
calculations is less than 1% of the total vectors, then 
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the data are considered reliable. 

Particle Seeding- Proper seeding of the flow 
is critical to accurate PIV measurements. The seed 
particles must be evenly distributed and of sufficient 
density to define the flow without altering its physical 
properties. The particles must also be small enough 
that they accurately follow the flow (especially chal- 
lenging in accelerating flows), yet large enough so that 
they scatter a sufficient amount of light to be de- 
tected. It is usually convenient to have remote control 
over the delivery of the seed material into the flow. 
The particle generator used in this test employed 
an inert gas to atomize a non-toxic, pharmaceutical- 
grade mineral oil. The mist was vaporized and then 
condensed before being released into the flow. The 
particle size was estimated to be less than 0.5 pm. 
The particle generator was located upstream of the 
settling chamber so that the particles would mix with 
the air entering the tunnel circuit before being drawn 
into the wake of the rotor. 

PIV Data Acquisition- The CCD cameras 
were connected to separate frame-grabber cards in- 
stalled in a PC workstation. Using Provision soft- 
ware, images were acquired and selectively processed 
immediately in order to evaluate the quality of the 
raw images and the adequacy of the inter-pulse time 
delay. Good image correlations depend on such is- 
sues as background light contamination, particle im- 
age brightness, contrast, focus, beam alignment, and 
light pulse separation. 

Prior to taking data, fine beam focusing and 
alignment were performed. Focusing was accom- 
plished by translating the camera lenses while ob- 
serving the real-time image displays. The gain and 
black levels were adjusted to maximize the contrast 
and brightness of the particle images. The laser sheets 
were judged to be coplan ar when both laser pulses 
were recorded as a single image and the resulting par- 
ticle images appeared as doublets. 

The correct inter-pulse delay can best be deter- 
mined under actual test conditions. Images were ac- 
quired and correlated to determine the maximum par- 
ticle displacement. This displacement corresponded 
to the maximum particle velocity and the pulse sepa- 
ration was adjusted to produce a ± 3 pixel range, re- 
gardless of the window size. Pixel displacements that 
are too small will limit the dynamic range whereas 
displacements that are too large will decrease the 
probability of correlation (ref. 7). Because of the high 
in-plane particle displacements, caused by the high 
circumferential velocity of the trailing vortex from the 
rotor at different wake ages, the optimum pulse sep- 
arations varied from 5 to 40 /is (tables 4 and 5). It 
should be noted that if the thickness of the light sheet 
is 1 mm and the tip of the rotor blade is moving at 


87 m/s, the time delay between laser pulses should 
be no more that about 12 /is to capture the high- 
est out-of-plane velocities just behind the rotor. If 
this condition is not satisfied at a calculation node, 
an interpolated value based on values at neighboring 
nodes results. Once image quality was assured, 25 im- 
age pairs were acquired per wake age. Data for each 
wake age required about 400 megabytes of storage. 

Table 2: Pulse Delay Times when 0 = 5° 


Wake Age 

Zone A 

Zone B 

Zone C 

0.50 c 

15 /zs 

20 jis 

12 /zs 

1.00 c 

15 /is 

20 /is 

20 /zs 

2.00 c 

15 /is 

20 /zs 

30 /zs 

3.00 c 

15 /is 

20 /zs 

40 /zs 

4.00 c 

15 /is 

20 /zs 

40 /zs 

5.00 c 

15 /zs 

20 ps 

40 /zs 


Table 3: Pulse Delay Times when 0=11° 


Wake Age 

Zone A 

Zone B 

Zone C 

0.50 c 

5 /zs 

5 /zs 

12 /zs 

1.00 c 

5 /zs 

5 /zs 

12 /zs 

2.00 c 

5 /zs 

5 /zs 

20 /zs 

3.00 c 

5 /zs 

5 /zs 

25 /zs 

4.00 c 

5 /zs 

5 /zs 

30 /zs 

5.00 c 

5 /zs 

5 /zs 

40 /zs 


Processing Wake Images- After completion of 
the test, the data were processed on a PC workstation 
using Provision software. To obtain good resolution 
of the primary vortex, a calculation mesh having 101 
nodes by 78 nodes was constructed over the region 
of significant interest in the flow field. The area cov- 
ered by the mesh resulted in approximately 17 pixels 
between nodes in both directions. The interrogation 
window was set at 24 pixels on a side, which gave 
a 70% overlap. The physical area measured in each 
image zone was 262 mm (with an axial increment of 
2.6 mm) by 149 mm (with a radial increment of 1.9 
mm). The 2.6 mm by 1.9 mm increments, which de- 
fined the calculation step sizes in the axial and radial 
directions, also defined the resolution of the measure- 
ments. 

A single file for each wake age includeed a flag 
for every velocity measurement that indicateed the 
nature of that vector (such as valid, invalid, inter- 
polated, recalculated, or not calculated). Reduced 
data files containing the coordinates of each calcula- 
tion node and the three components of velocity (in 
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terms of displacement) were stored in ASCII format 
and then transferred to a mainframe computer for 
analysis. 

The velocity components were first converted 
from displacement units to velocity units based on 
the pulse duration for that particular measurement. 
The physical coordinates were transformed so that 
r = z = 0 would correspond to the center of the rotor 
hub (fig. 7). The sequence for extracting the char- 
acteristics of the trailing vortex (fig. 8) for each case 
began with a calculation of the vorticity field, for 
every image pair. Vorticity was calculated by integra- 
tion (circulation box method) according to: 

u>e = A“ 1 F s= A” 1 <j> (u, t?) • d(r } z) (1) 

where the direction of integration was such that the 
enclosed area, A, was on the left of the integration 
path. The size of each circulation box was defined 
by the local mesh spacing, which on average rendered 
A = 5 mm 2 . 

When more than one vorticity center was present 
in the field of view, a single image pair had to be pro- 
cessed to determine where the field should be parti- 
tioned so that statistics could be built on only one 
vortex. A simple circular partition that focused on 
the vortex of interest was usually sufficient to sepa- 
rate the two vortices. 

After calculating the vorticity distribution over 
the entire image, a search was performed to discover 
the zone of greatest vorticity concentration based on 
values exceeding 50% of the maximum. The center of 
vorticity for a given image was defined to be at the av- 
erage of these selected locations. The 50% threshold 
criterion prevented the center of vorticity for an im- 
age from being influenced by a single, and possibly er- 
roneous, measurement of velocity (which in turn con- 
tributes to the calculation of vorticity). The resulting 
centers of vorticity (one for each image) formed a set 
for which the mean value and the standard deviation 
could be calculated. Normally any image with a vor- 
tex center that deviated from the mean by more than 
1.5 times the standard deviation would be purged 
from the set. Since the flow in the near wake was ex- 
ceptionally periodic in this experiment, the resulting 
standard deviation was small enough that all images 
were eligible for inclusion in the average. The stan- 
dard deviation criterion prevents an atypical depar- 
ture, albeit a legitimate image, from being factored 
into the average. The threshold and purging criteria 
both have the affect of reducing the degree of wander 
in the data. 

The images must be averaged to smooth out 
small irregularities in the flow. However, whenever 


there are features in the flow with high gradients, such 
as a trailing vortex, and the vortex appears in differ- 
ent positions from image to image (vortex wander), 
important details of the structure tend to be smeared 
out if a simple average of the images is performed. 
The remedy was to perform the average after artifi- 
cially aligning the images based on a recurring feature 
of the flow, which in this case was the center of vor- 
ticity. This process, in contrast to a simple average, 
is referred to as conditional ensemble averaging; the 
condition here being the alignment of all the images in 
the set based on their centers of vorticity. The image 
with a vortex center closest to the mean location that 
was established for the set of images was selected as 
the “anchor” . The indices of each image matrix were 
then adjusted according to the offset of each vortex 
center from the vortex center in the anchor image. 
All data in the anchor image were retained. However, 
some portion of all other images with adjusted indices 
that fall outside the boundaries of the anchored image 
were necessarily discarded. Hence, as a consequence 
of vortex wander, the population contributing to the 
conditional average was greatest over the interior of 
the matrix. 

To extract information about the geometry of the 
vortex (see ref. 8 for more detail), the center of the 
vortex was assumed to be at the center of swirl (which 
was not necessarily at the center of the vorticity that 
was accreted from the rotor wake, which had the form 
of a spiraling sheet with varying vorticity). Using the 
center of vorticity as a starting point, the surround- 
ing locations were interrogated to determine the best 
node for which the sum of the dot products of two 
unit vectors was a minimum over a neighborhood of 
locations surrounding the candidate node. The coor- 
dinates of the neighboring node relative to the candi- 
date node defined one of the unit vectors. The veloc- 
ity vector at the neighboring node defined the other 
unit vector. Denoting the angles of these two unit 
vectors by a and /? (fig. 9), then the node nearest the 
center of swirl was found when the following expres- 
sion was a minimum: 

N 

Y2 |cos(or- 0)\/N (2) 

n= 1 

where N is the total number of neighboring points 
considered. Relative to the center of swirl for a pure 
vortex in a stationary flow, a and /? will be orthogonal 
for all neighboring points and the sum will be exactly 
zero. For most real flows the sum will not be precisely 
zero, and the procedure works best when the neigh- 
borhood under consideration does not extend beyond 
about one or two core diameters (which is normally 
sufficient to cover a majority of the vorticity in a sin- 
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gle trailing vortex) 

After locating the center of swirl, the size of the 
vortex core can be found by dividing the surround- 
ing region into annular zones and calculating the av- 
erage swirl velocity and the associated mean radius 
of each annulus. This data can then be fit with a 
least-squares spline subject to the condition that the 
resulting curve have only one inflection along its in- 
ner extent (nominally set to twice the estimated core 
radius that was previously found). The point of inflec- 
tion now defines the radius of the vortex core, since 
that is where a maximum value for the swirl veloc- 
ity is reached. This procedure was considered to offer 
a more rational approach for determining the size of 
the vortex, especially when the velocity peaks that 
are characteristic of vortices (and upon which core 
sizes are traditionally based) are dependent on how 
the vortex is sliced (vertical, horizontal, or otherwise). 

Test Conditions- The experiment was per- 
formed at a constant rotor speed of 870 RPM (14.5 
Hz), which for a rotor diameter of 75.4 inches cor- 
responds to a tip speed of 286.2 ft/s. Based on an 
average ambient temperature of 62° F, the Reynolds 
number based on tip chord was Re = 0.6 x 10 6 and 
the Mach number was M = 0.26. The collective pitch 
angle and the rate of climb were the only variables in 
this test (table 1). The lowest climb condition was in- 
tended to occur when the return flow was not blocked 
and the fan was idle. This is referred to as the “nor- 
mal tunnel” configuration, with the rotor itself sub- 
stituting for the tunnel fan to establish a modest flow 
around the tunnel circuit. 

As indicated in table 1 , wake measurements were 
performed only for collective pitch angles of 5° and 
11°. In order to obtain high resolution PIV images, 
the cameras were focused on a small region in the near 
wake behind the trailing edge of the rotor. To expand 
the physical area under study the camera view was ra- 
dially translated over a total of 3 contiguous regions, 
giving a total coverage of the outer 56% of the ro- 
tor blade. Measurements at different wake ages were 
obtained by adjusting the timing of the pulsed light 
sheets, thereby permitting the tip of the rotor blade 
to incrementally pass beyond the plane of the light 
sheets by up to 5 chord lengths. The PIV measure- 
ment locations are shown in table 2. 

RESULTS AND DISCUSSION 

Equilibrium Conditions- Given that it takes a 
certain amount of time for the flow in the tunnel to ad- 
just to the action of the rotor and to reach equilibrium 
(during which time balance zeros may drift due to 
heating), time histories of the rotor loads and the tun- 
nel velocity were recorded (fig. 10). The results indi- 
cate that the velocity in the test section and the rotor 


Table 4: Test Conditions 


Collective Pitch Angle 


Configuration 

3° 

5° 

7° 

9° 

11° 

Tunnel Blocked 

y/ 

* 

yj 

yj 

* 

Normal Tunnel 

V 

y/ 

V 

V 

V 

Climb, fan @ 22 RPM 

V 

V 

V 

V 

V 

Climb, fan @ 52 RPM 

V 

V 

y/ 

V 

V 

Climb, fan @ 80 RPM 

V 

V 

yf 

V 

V 

Climb, fan @ 105 RPM 

V 

V 

V 

V 

V 


* including PIV wake measurements 


Table 5: Location of PIV Measurements 


Wake Age 

Zone A 

Zone B 

Zone C 

0.50 c 

V 

V 

y/ 

1.00 c 

V 

yj 

yj 

2.00 c 

V 

yj 

yj 

3.00 c 

V 

yj 

yj 

4.00 c 

V 

■J 

yj 

5.00 c 

V 

yj 

yj 


torque reached acceptably level values within 5 min- 
utes. However, the thrust (which was measured with 
a 6-component balance) asymptotically approached a 
steady value over a period of about 20 minutes. This 
suggests that the air in the tunnel circuit was set into 
motion rather quickly and that the continued change 
in thrust was due to balance heating. 

Since the balance appears to be considerably 
more sensitive to temperature at low thrust levels, it 
is clear that a pre-run tare reading will not be useful. 
To evaluate the suitability of a post-run zero, time 
histories of the thrust and torque were recorded dur- 
ing successive start ups and shut downs of the rotor 
(fig. 11). Using arbitrary references, the buildup and 
decay in thrust readings are distinctly different from 
the torque readings. After a shutdown was initiated, 
the rotor quickly came to rest while the flow in the set- 
tling chamber (where the rotor was located) became 
essentially zero after one minute. During that brief 
one-minute time interval the diminishing impact of 
the air stagnating on the blades appears to be some- 
what offset by the temperature drift of the balance. 
Although not ideal, the thrust and torque readings 
at one minute after shutdown were taken to be the 
correct tare values. 

Rotor Performance- The wake of the rotor was 
free to locally recirculate and become reingested by 
the rotor, while the remainder of the wake passed 
downstream and through the test section. As de- 
picted in the upper contour plot (fig. 12), the flow 
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in the test section was not only nonuniform, but fa- 
vored one side. The skewing of the flow toward one 
side is believed to be due to the blockage effect of the 
PIV support structure on the opposite side of the tun- 
nel, which restricted the amount of air flowing near 
the walls of the settling chamber as it was being en- 
trained by the rotor wake and drawn into the test 
section. Adding an additional 37% to the area sur- 
veyed (lower contour plot) yielded only a negligible 
change in the averaged velocity. 

With the return flow blocked, rotor loads were 
measured for collective angles 3° to 11° (fig. 13). The 
subplot showing the variation in test section velocity 
with collective pitch angle shows that the single mea- 
surement at the centerline of the tunnel reasonably 
approximates the value obtained by averaging over 
the central portion of the test section. Using these 
values to represent the flow through the test section, 
conservation of mass and the known area ratio dictate 
the velocity that must exist upstream of the rotor. 
These deduced values agree well with the velocities 
that were measured upstream of the rotor (about 4.5 
rotor diameters). Based on thrust and torque mea- 
surements, the calculated hover figure-of-merit (FM) 
is also shown for the range of collective pitch angles. 

FM = {T*f2pA) l f*/Qn (3) 

It should be noted that the above relation assumes 
that the flow is unbounded and that the far upstream 
velocity is zero. Although these assumptions were vi- 
olated in this experimental setup, especially when the 
tunnel fan was used to create a climb condition, the 
hover figure-of-merit will simply serve as a convenient 
nondimensional group for comparison purposes. 

With the tunnel circuit in its normal configura- 
tion (return flow unblocked), that portion of the ro- 
tor wake that is not skimmed off near the walls of 
the heat exchanger section will continue around the 
tunnel. That means that the rotor had to be drawing 
stagnant air from outside the facility and mixing it 
with a portion of the wake that the rotor had already 
generated. In essence, even with the tunnel fan at 0 
RPM, the rotor was in a slight climb condition. With 
all other parameters remaining the same, the result- 
ing rotor loads were slightly lower (fig. 14). The mea- 
sured velocity approaching the rotor was perceptively 
higher, a condition that was made more evident in the 
“amplified” flow that developed in the test section. 

Running the tunnel fan at speeds from 22 to 105 
RPM increased the upstream velocity until it com- 
pletely dominated the contribution of the rotor to the 
recirculating flow in the tunnel circuit. The load vari- 
ation with collective pitch angle retained the same 
shape (figs. 15-18), however the loads decreased with 


increasing “climb rate” due to an effective reduction 
in the blade incidence (fig. 19). The climb velocity 
is defined as the relative motion of the rotor disk to 
free air. Since the tunnel walls bound the rotor, the 
upstream velocity was not entirely the climb velocity. 
However, for simplicity, the climb rate in this paper 
is taken to mean Voo/Vup , where Voo refers to the 
upstream velocity and V t i p = RQ. 

It would appear that a good estimate for the 
hover figure of merit can be obtained by extrapolat- 
ing the curve for any given collective pitch back to 
Voo /Vtip = 0 (as suggested in ref. 2). Although the 
placement of the rotor in the settling chamber of the 
wind tunnel did result in a smoother flow than would 
normally be achieved in the Hover Test Facility at 
Ames Research Center (a comparatively small cubic 
space with an asymmetric flow intake), the confin- 
ing boundaries of the wind tunnel do affect the flow 
and may prevent the simulation of a true hover con- 
dition, even with the return flow blocked. When the 
return flow in the tunnel is blocked and the rotor 
wake that passes through the test section is totally 
expelled, continuity requires that a flow upstream of 
the rotor must exist, and this contradicts the ideal 
hover state. Given that an upstream flow did exist 
in this experiment, it is difficult to imagine how this 
flow could have passed through the rotor entirely, and 
not have at least a portion of the flow travel around 
the rotor. By considering the case when the return 
flow was blocked, one can show that a portion of the 
flow did indeed bypass the rotor. If the rotor were to 
have been operating in free air, it can be shown from 
momentum theory that: 



where V a represents the flow through the rotor disk, 
and A represents the area of the disk (31 ft 2 ). For 
9 = 11° the thrust is about 40 lb/. The volume flow 
rate through the rotor disk is AV a = 505 ft 3 /s. How- 
ever, based on the averaged velocity measured across 
the test section, the volume flow rate passing through 
the tunnel is 1444 ft 3 /s, which is about 3 times higher 
than the volume passing through the disk. This im- 
plies that the rotor had set into motion a great deal 
more fluid upstream of the rotor than actually passed 
through the rotor disk, and places some degree of un- 
certainty on just what flow state can be ascribed to 
this experimental setup. 

As the wake of the rotor (along with the sur- 
rounding air that had been entrained) passed through 
the tunnel contraction, the entire flow was accelerated 
and the static pressure decreased. This difference 
in pressure between the test section and the settling 
chamber must have lead to the development of an up- 
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stream velocity that appeared to the rotor as a climb 
velocity. This unintentional climb velocity could have 
been forced to zero by restricting (or blocking) the 
flow through the test section, but then the rotor would 
have been operating in a ” ground effect” state. A 
better solution may be to eliminate this pressure dif- 
ferential by exposing the test section to atmospheric 
pressure, although the tunnel contraction boundaries 
would undoubtedly still have some influence on the 
flow. A numerical analysis was performed on a rotor 
in free air as well as with straight wall and tunnel 
contraction boundaries (ref. 8). These results indi- 
cate that the tunnel walls do indeed have an effect 
on the flow and that the correct hover figure of merit 
(which was shown to be higher) can only be obtained 
in free air. 

Rotor Wake- The following discussion of the 
rotor wake will initially focus on the tip region of 
the flow (zone A) over progressively larger wake ages 
(from ip = 3° to 31°). Each series, first for 0 = 5° 
and then for 0 = 11°, will consist of two parts: (a) 
the location and size of the vortex and (b) contour 
maps of the velocity and vorticity fields. 

Statistics on the location and size of the vortex 
(part (a) of each figure) includes a graphic showing 
the degree of wander for the set of images (for exam- 
ple, fig. 20). The grid spacing of A r — 1.9 mm and 
Az = 2.6 mm corresponds to the calculation mesh 
used to determine the velocity field. The vorticity 
field is shown as a color contour map. The magnitude 
of the vorticity is represented by the level of color sat- 
uration, with a blue hue denoting u/$ < 0 and a red 
hue denoting > 0. Along with the position of the 
rotor tip, the location of the vorticity center (defined 
at the centroid of vorticity values exceeding 50% of 
the maximum) is also shown. The inplane velocity is 
shown as colored vectors with the hue ranging from 
blue (zero velocity) to red (maximum velocity). The 
location of the center of swirl (as defined earlier) is 
also shown. The circled area identifies the range of 
swirl vectors used to construct the swirl velocity fit, 
which in turn defines the peak swirl velocity and the 
size of the vortex core. 

The velocity and vorticity contours (part b of 
each figure), which include a superposition of contour 
lines, more clearly emphasize the relationship between 
the tip vortices and the remainder of the rotor wake. 
Only blue and red hues are used, with the magnitude 
of the variable indicated by the saturation of the two 
colors. Maximum and minimum values of each vari- 
able are shown along side of the respective color bars. 

The results for 0 = 5° (figs. 20 - 25) show that 
the vortex wandered by an average of 0.65 mm, which 
is about 8% of the core radius. The vortex from the 
second blade (ip = 180° earlier) accompanied the vor- 


tex from the first blade in all zone- A images. The in- 
board portion of the wake that is associated with the 
first vortex appears to be distorted by the close prox- 
imity of the second vortex at all wake ages studied (up 
to ip = 31°). Aside from the tip vortex, the wake from 
the second blade cannot be detected. The difference 
in the locations of the vorticity center and the swirl 
center averages about 2 to 3 mm, or about 32% of the 
core radius. The viscous wake along the span of the 
first rotor blade is clearly evident in the strong out-of- 
plane velocity component (u), which shows a notice- 
able disconnect in magnitude just as it approaches the 
tip vortex. 

The results for 0 = 11° (figs. 26-31) show that 
the vortex wandered by an average of 0.54 mm, which 
is less than in the 0 = 5° case, yet is about 9% of 
the core radius that occurred at this higher collec- 
tive angle. The axial separation between the vortices 
from the first and second blades has substantially in- 
creased. Although the axial position of the viscous 
wake from the first blade has also increased, the dis- 
tance between this wake and the second vortex is now 
sufficient (compared to the 0 = 5° case) that no dis- 
tortion of the wake appears to be present. Aside from 
the tip vortex, the wake from the second blade can- 
not be detected (similar to the 0 = 5° case). The 
difference in the locations of the vorticity center and 
the swirl center averages little more than 1 mm, or 
about 21% of the core radius. The out-of- plane ve- 
locity component continues to show a disconnect be- 
tween the flow in the tip vortex and the remainder of 
the viscous wake. 

The data does not support any obvious trend in 
either the size or intensity of the vortex during the 
first 31° of wake age (fig. 32). On average, r c = 7.7 
mm (7.4% of c) when 0 = 5° and decreases to r c = 6.3 
mm (6.1% of c) when 0 = 11°. The maximum swirl 
velocity averages 7.4 m/s (8.5% of RQ) when 0 = 5° 
and increases to 29 m/s (33% of RQ) when 0 = 11°. 
The radial location of the tip vortex decreases rapidly 
with wake age, traveling in a linear fashion at about 
0.5 mm per degree of wake age (totaling 1.5% of R) 
when 0 = 5° and nearly double this amount when 
6 = 11°. The axial location of the tip vortex is more 
difficult to summarize. In the 0 = 5° case, the tip 
vortex appears to suddenly jump in elevation after 
ip = 10° and remain at this level. On the other hand, 
in the 0=11° case, the tip vortex travels in an in- 
creasingly steep axial descent with a total displace- 
ment of about 1.1% of R after reaching tp = 31°. 

The extent of the imprint along the wake of the 
rotor blade can be more fully appreciated by assem- 
bling all 3 measurement zones. In the 0 = 5° case 
(fig. 33), this presentation encompasses the position 
of the tip vortex after 1 revolution (zone B). Clearly 
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the bulk of the vorticity in the wake is primarily con- 
fined to the tip vortices. It is also evident that a 
deficit in the out-of-plane velocity component (t;) ex- 
ists all along the wake of the blade as well as within 
the tip vortices. At this wake age (tp = 12.4°) the 
maximum magnitude of the out-of-plane component 
is about 54% of the maximum swirl velocity. 

Although dwarfed by the intensity of the tip vor- 
tex, selected regional enlargements and a rescaling of 
the color contours can enhance the appearance of vor- 
ticity that is present in the remainder of the wake. In 
the 0 = 5° case (fig. 34), a thin stretch of positive 
vorticity overlays an equally thin stretch of negative 
vorticity over a majority of the wake. The positive 
vorticity adjacent to the first vortex (fig. 34a) proba- 
bly relates to the change in blade circulation and is in 
the process of being accreted into the primary vortex. 
The patches of negative vorticity located above and 
to the left of both the first and second vortices are due 
to the shear layers created between these vortices and 
the surrounding fluid (fig. 35). The layered pattern of 
positive and negative vorticity along the inner portion 
of the blade wake (figs. 34b and 34c) may relate to the 
shear layers caused by a centrifugal pumping of the 
fluid within fche boundary layers of the rotating blade 
(thereby causing a radially outward flow on both the 
upper and lower surface of the blade). Another reason 
may be the changing strength of the spanwise vortic- 
ity in the upper and lower surface boundary layers, 
which would in turn require the shedding of stream- 
wise oreinted vorticity (positive and negative). 

Similar flow patterns exist in the 5 = 11° case 
(figs. 36 and 37). However, one distinctive difference 
is the greater separation between the first and second 
tip vortex. As a result, the blade wake associated 
with the first vortex does not appear to be distorted 
by the second vortex. Although the second vortex 
continues to have some affect, theories attempting to 
relate measured wake vorticity to conditions on the 
blade can be expected to be more accurate in cases 
where the wake is convected away quickly and the in- 
teraction between the wakes from different blades is 
minimized. Another difference between the two col- 
lective pitch cases is the maximum magnitude of the 
out-of-plane component, which for the 9 = 11° case 
is only about 24% of the maximum swirl velocity. 


Table 6: Vortex Characteristics at xp = 12.4° 


e 

r c /c 

(it + v)/R£t 

w/Rn 

5° 

6.4% 

8.9% 

8.4% 

n° 

6.1% 

33.8% 

11.0% 


Vortex characteristics at a distance of 2 chord 
lengths into the wake are presented in table 6. Nondi- 
mensional groups are formed using the chord, radius 
and rotational speed at the tip of the rotor blade. 

CONCLUSIONS 

1. Although the confining boundaries of the 
facility prevented ideal hover conditions from being 
achieved, a steady flow environment was formed in 
which rotor performance and wake structures could 
be studied at low rates of climb. 

2. The hover figure of merit increased with in- 
creasing collective pitch angles. A reasonably well- 
behaved and steady increase in this quantity was ob- 
served as the climb rate diminished toward zero. 

3. A majority of the vorticity in the wake is con- 
centrated in the tip vortices. The absence of any size- 
able levels of vorticity (within an order of magnitude) 
was noted in the remainder of the wake, even at the 
early wake age of \p = 12°. 

4. Aside from the persistent presence of the tip 
vortices generated by each passing blade, the imprint 
of the blade on the remainder of the wake is evident 
primarily in terms of the out-of-plane components of 
velocity and vorticity. However, this evidence appears 
to fade quickly and cannot be observed at all in the 
wake of the following blade (ip — 180°). 

5. When the rotor thrust is low (for 9 = 5°), 
the inboard portion of the wake is distorted by the 
close proximity of the tip vortex from the preceding 
blade. The axial separation distance is greatly ex- 
tended when the rotor thrust is increased (0 = 11°), 
which may allow the vortical composition of the near 
wake to be more easily identified with the drag com- 
ponents (induced and profile) on the rotor. 

6. The size of the vortex core in the near wake 
(rp < 31°) is about 6% of the blade chord for col- 
lective pitch angles of 5° and 11°, At ip = 12°, the 
out-of-plane velocity component is about 8% of the ro- 
tational speed of the blade tip at 9 = 5° and the swirl 
velocity is only slightly higher. However, at 9 = 11°, 
the out-of-plane velocity increases to 11% of the ro- 
tational speed of the blade tip, but the swirl velocity 
nearly triples the magnitude of the out-of- plane ve- 
locity. 
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Figure 1: 7- by 10-Foot Wind Tunnel flow circuit. 
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Figure 2: Instrument setup in settling chamber of 7- by 10-Foot Wind Tunnel. 
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Figure 3: OLS blade profile. 
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Figure 4: OLS Model rotor. 
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Figure 6: Perspective effect of particle displacement with Stereoscopic PIV. 
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Figure 7: Particle imaging system layout. 
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Figure 8: Sequence for extracting vortex characteristics. 
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Figure 9: Angles used in search for center of swirl. 
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Figure 10: Time required to reach equilibrium conditions. 
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Configuration: Tunnel Return Flow Blocked 
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Figure 12: Velocity distributiou in test section with return flow blocked. 
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Air return section blocked ... fan speed = 0 RPM 
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Test section velocity ... view looking downstream 





Figure 13: Rotor loads and tunnel velocities with return flow blocked. 














Normal wind tunnel configuration ... fan speed = 80 RPM 
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Test section velocity ... view looking downstream 
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Figure 18: Rotor loads and tunnel velocities with return flow unblocked and fan speed = 105 RPM. 
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Figure 20: Near wake at 9 = 5° and xj) = 3.1° (0.5 chord downstream from rotor tip). 
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Figure 20: Concluded. 
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(a) Location and size of vortex. 

Figure 21: Near wake at 0 = 5° and ip = 6.2® (1 chord downstream from rotor tip). 
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Figure 22: Near wake at 9 = 5° and ip = 12.4° (2 chords downstream from rotor tip). 
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Figure 23: Near wake at 6 — 5° and tp = 18.6° (3 chords downstream from rotor tip). 
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Figure 23: Concluded. 
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(a) Location and size of vortex. 

Figure 24: Near wake at 6 = 5° and ip = 24.8° (4 chords downstream from rotor tip). 


38 








Swirl velocity, m/s 


Ar= 1.9 mm 


100 


Vortex 

wander 


Az= 2.7 mm 


50 h 


0 

E 

E -50 

N 

-100 


Rotor Tip 


-160 h 


Collective*^ 0 Wake Age“5.0 chords 
std dev allowed= 3.5 Zone: A 
account for wander: focus 
max u>q threshold= 60 % 


935 mm 


1 Vorticity 

Center 

( -2 mm 


I 


! — ► I 

\<*B\ 


1/s 

5487 



-5487 



r # mm 



Max swirl* 3 7.6 m/s 
Core radius= 7.3 mm 


(a) Location and size of vortex. 

Figure 25: Near wake at 6 = 5° and ip = 31.0° (5 chords downstream from rotor tip). 
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Figure 26: Near wake at 9 = 11° and ip = 3.1° (0.5 chord downstream from rotor tip). 
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Figure 27: Near wake at 0 = 11° and ip = 6.2° (1 chord downstream from rotor tip). 
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(b) Velocity and vorticity contour maps. 
Figure 27: Concluded. 
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(a) Location and size of vortex. 

Figure 28: Near wake at 6 = 11° and rp = 12.4° (2 chords downstream from rotor tip). 
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(b) Velocity and vorticity contour maps 
Figure 28: Concluded. 
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Figure 29: Near wake at 8 = 11° and rp = 18.6° (3 chords downstream from rotor tip). 
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(b) Velocity and vorticity contour maps 
Figure 29: Concluded. 
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(a) Location and size of vortex. 

Figure 30: Near wake at 0 = 11° and ip — 24.8° (4 chords downstream from rotor tip). 
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(b) Velocity and vorticity contour maps 
Figure 30: Concluded. 





Collect jve=11° Waka Age=5.0 chords 
std dev allowed^ 3.5 Zone: A 
account for wander: focus 
max uq threshold^ 50 % 


923 mm 



Vorticity 

Center 

- - -9 mm 




925 mm 



500 600 700 800 900 1000 

r, mm 



o 10 20 30 40 60 

r, mm 


(a) Location and size of vortex. 

Figure 31: Near wake at 6 = ll® and ip = 31.0° (5 chords downstream from rotor tip). 
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(b) Velocity and vorticity contour maps. 
Figure 31: Concluded. 
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Figure 33: Velocity and vorticity contours covering 3 radial zones of flow for 0 = 5° and if) = 12.4° (2 chords) 
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Figure 35: Enlargement of inplane velocity in zone A for 9 = 5° and ^ = 12.4° (2 chords). 
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